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ABSTRACT 

HI features near young star clusters in M81 are identified as the photodissociated surfaces 
of Giant Molecular Clouds (GMCs) from which the young stars have recently formed. The HI 
column densities of these features show a weak trend, from undetectable values inside R = 3.7 
kpc and increasing rapidly to values around 3 x 10^^ cm~^ near R k, 7.5 kpc. This trend is 
similar to that of the radially-averaged HI distribution in this galaxy, and implies a constant area 
covering factor of k, 0.21 for GMCs throughout M81. The incident UV fluxes Go of our sample 
of candidate PDRs decrease radially. 

A simple equilibrium model of the photodissociation-reformation process connects the ob- 
served values of the incident UV flux, the HI column density, and the relative dust content, 
permitting an independent estimate to be made of the total gas density in the CMC. Within the 
CMC this gas will be predominantly molecular hydrogen. Volume densities of 1 < n < 200 cm~^ 
are derived, with a geometric mean of 17 cm~^. These values are similar to the densities of GMCs 
in the Galaxy, but somewhat lower than those found earlier for MlOl with similar methods. Low 
values of molecular density in the GMCs of M81 will result in low levels of collisional excitation 
of the CO(l-O) transition, and are consistent with the very low surface brightness of CO(l-O) 
emission observed in the disk of M81. 

Subject headings: galaxies: individual: M81 - galaxies: ISM - ISM: clouds - ISM: H2 - galaxies: FUV - 
galaxies: HI 



1. Introduction 

Giant molecular clouds (GMCs) in the inter- 
stellar medium (ISM) are generally accepted as 
the birthplaces of new stars. The most mas- 
sive of these new stars produce copious amounts 
of far ultraviolet (FUV) radiation which will, 
in turn, photodissociate the molecular gas of 
the parent GMCs, producing "blankets" of HI 
(and other atomic specie s ) on their surfaces. 
Allen. Atherton fc Tilanui (|l986f) were the first 
to present evidence that major features in the 
HI distribution of the nearby spiral galaxy M83 
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(NGC 5236), namely, the inner HI spiral arms 
observed with the Very Large Array (VLA), were 
the result of photodis s ociati on of H2 on galac- 
tic scales. Allen et al] ()l997t ) confirmed that HI 
features existed in M81 (NGC 3031) on scales of 
w 150 pc which were qualitatively consistent with 
the expected morphology of large, low density 
photodissociation regions (PDRs), and exphcitly 
related those HI features to nearby bright sources 
of Far-UV (FUV) radiation found on images of 
the galaxy from the Ultraviolet Imaging Telescope 
( UIT). 

ISmith et~aL ( 2000l ) used a simple, but quanti- 
tative, model for the equilibrium physics of pho- 
todissociation regions and applied it to VLA-HI 
and UIT-FUV observations of MlOl (NGC 5457). 
Their work showed that this approach provides an 
entirely new method for determining the volume 
density of molecular gas in star-forming GMCs of 
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galaxies, a method which is no longer dependent 
on the (often poorly known) excitation conditions 
for line emission by specific molecular tracers. 

In this paper we return to another spiral galaxy, 
M81, and carry out a quantitative analysis of the 
G MCs in this galaxy using the methods described 
bv lSmith et"al](l2000l) . This new analysis has been 
made possible by new data on M81 which was 



not available to Allen et al. ([1993); 

new, higher- 



resolution, high-sensitivity VLA-HI observations, 
and sensitive new FUV imagery from the GALEX 
satellite. We have also sought an independent 
verification that the very basis of our approach 
is valid, namely, that the HI in the immediate 
vicinity of FUV concentrations is indeed produced 
in PDRs. To this end we have examined the 
Spitzer/IRAC data on M81 for evidence of mid- 
IR emission by polycyclic aromatic hydrocarbons 
(PAHs) which are also thought to be tracers of 
PDRs. 

M81 has been probed extensively for the molec- 
ular tracer CO. No global CO map of M81 has 
been published to this date, and the detected CO 
emission is found to be very w eak and spotty 



Knapen et al. (l2006h and references 



(see e.g. 

ther ein). Since PDRs also produce CO emission 
(e.g. lAllen et al.ll2004r ). comparing the CO results 
to the results in this paper is therefore of consid- 
erable interest. 

The outline of this paper is as follows: Section 
[2] contains a description of the data we used, fol- 
lowed by a brief theoretical description of PDRs 
in Section [3] and the application of the method in 
Section m The results are presented in Section [H 
The results are discussed and the conclusions are 
briefly summarized in Section [S] 

2. Data 

Th e distance to M 8 1 is ta ken to be 3.6 Mpc, 



after Freedman et al 



tion of 



58° 



(de Vaucouleurs et al 



with an inclina- 
and positio n angle P. A. ~ 157° 
9911 ). At this distance, 
plane of the sky. The 
in all the plots have been 



1" = 17.5 pc in the 
galactocentric radii 
normalized to i?2R, which is 1 2.66 kpc for M81 
(|Vila-Costas fc Edmund^ Il992l ). For the model 



we used (Section [3]), we need HI data, FUV data 
and the dust-to-gas ratio scaled to the local solar 
neighborhood value {6/6q). 



Gamett 5, Shields (1987) 



Stauffer S. Bothun (1984) 



Fig. 1. — Dust-to-gas ratios as a function of R 
in M81, as derived from lCarnett fc Shield^ (|l987l ) 



(dashed line) and lStauffer fc BothunI (|l984l ) (solid 
line). So is the value for the solar neighborhood in 
the Galaxy. We use the Stauffer & Bothun data 
in this paper. 

A new, high quality integrated HI column den- 
sity map of M81 made with the V LA has been 



made available by the THINGS team ([Walter et al. 
20051 ). It has an angular resolution of 7.58" x 



7.45", equivalent to about a 130 pc linear resolu- 
tion. The pixel spacing is 1.5". An azimuthally- 
averaged radial profile of this data has been made 
available by Frank Bigicl (THINGS, private com- 
munication) . 

The ultraviolet data for M81 has been obtained 
from the publicly accessible GALEX data archive. 
The angular resolution is 4" at a pixel spacing of 
1.5". The astrometry has been checked against 
a number of stars in the Guide Star Catalogue 
v2.3. The Near-UV image was used for this pur- 
pose, since not all of these stars are detected in the 
Far-UV image. The NUV and FUV images were 
verified to have the same astrometry. We use the 
FUV image at A « 150 nm as a "proxy" measure 
of the photons with A < 100 nm that can directly 
dissociate molecula r hydrogen, as was done in e.g. 



Smith et all (|2000h 



The dust-to-gas ratio {S/Sq) is derived from the 
metallicity measure 12 -1- log{0/H). In the case of 
M81, the slope of log d /6n is taken as equal to the 



slope of 12-flog(0/iif ) (jlssa. MacLaren fc Wolfendale 
1990l ). In the solar neighborhood. 12+ l og(0/ g) = 
8.8, according to iMihalas fc BinnevI (|l98lh . At 
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FUV HI IR 1 




R.A. (2000) 



Fig. 3. — Detail plots of four sources from Table [T] Their galactocentric radius is 4.1, 6.1, 7.8 and 9.6 kpc 
respectively. All sources are included in the electronic version. The HI is shown as the underlying grayscale 
in every panel. The central panel also shows HI contours. The lowest and highest contours are colored black 
and white respectively. The UV source (star) is marked in every panel. Left panel: UV contours (purple) 
start at 6 X 10~^® ergs cm~^ s~^ with each consecutive contour 1.48 times the previous one. The 

HI grayscale is as in the middle panel, but at the original pixel size of 1.5". Middle panel: HI grayscale 
following the contours, of which the starting contour level and increment (previous level times increment) 
are given in Table[T] Right panel: Non-stellar (PAH, red) contours at relative levels, starting at 1.25%, with 
each consecutive contour 1.58 times the previous one. The HI grayscale is as in the middle panel, but at 
the original p ixel size. Note th e HI surrounding the UV sources in 'blankets', which becomes clearer after 
zooming out (jAllen et al.lll997l) . 
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R.A. (2000) 

FUV HI IR 21 




R.A. (2000) 

Figure [3]- continued. 
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Fig. 4. — HI projection issues. Panel a: MSl's position angle and inclination are used to deproject pni- 
Panel b: the exact vertical position of the HI patch along the line of sight is unknown. Panel c: the measured 
HI column is assumed to be a sufficient measure of the HI column needed for the PDR model, although the 
geometry of the HI layer is unknown. 



this value, S/Sq equals 1 by definition. This is suf- 
ficient to derive the expression we need for S/Sq- 

The 12 + l og(0/g) data for M81 has been 
adopted from IStauffer fc BothunI (Il984l). but we 
also c onsidered the results of Garnett fc Shield^ 
(jl987l ). The gradient they found is likely to be on 
the high side (Garnett, private communication). 
Figure [T] shows both 6 /So fits and illustrates the 
possible spread in the gradient oi S/Sq, similar to 
the spread i n 12 + \or(0/H) t h at is implied from 
Fig ure 4 of iGarnett fc Shieldsl (|l987l ). The fit to 
the Stauffer fc Bothun (jl984l ) data we used yields 



log(<5/^o) 



-0.045i?+ 0.198, 



(1) 



where R is the relevant galactocentric radius in 
kpc for each region under consideration. 

In order to identify PDRs from their PAH emis- 
sion, we used M81 S pitzer data from the SINGS 
public release (Kennicut et al.l 120031 ). In particu- 



lar, we used their 3.6 /um and 8.0 ^m data from 
the IRAC camera. To verify that this data lined 
up correctly with our other data, we checked the 
astrometry against the same GSC stars as we used 
for the UV data and found the astrometry to be 
accurate to within 1". 



3. Theory 

We use the data described in the previous sec- 
tion to derive the desired physical quantity n, the 
total density of gas in the GMC; n — nni + 
2nff2 (cm"'^). The atomic hydrogen column den- 
sity, FUV flux and dust-to-gas ratio are the re- 
quired parameters to calculate the total hydrogen 
volume density. 

Atomic hydrogen is slowly but continuously 
converted to molecular hydrogen in the interstel- 
lar medium, a process that is catalyzed by dust 
grains. Molecular hydrogen is rapidly photodisso- 
ciated back into atomic hydrogen by UV photons. 
Simple models for equilibrium concentrations of 
these component s have b een discuss ed by sev- 



eral a uthor s, e.g. ISternberg ( 1988.) , Wolfire et al 
(|l990l) and lHollenbach et all |l99l . 



Alien! (120041 ) has reviewed the physics and com- 



pared the time scale for equilibrium with other 
time scales i n galactic disks. His Equation 6 as 
derived from ISternbergl ( 1988h explicitly includes 
the dependence on the dust-to-gas ratio, as well as 
im proved values of th e coefficients as determined 
by I Allen et"al] (|2004[ ) in t heir Equation A2. The 
equation we shall use from I Allen (2004) is: 
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Fig. 5. — Individual (not ring-averaged) HI col- 
umn densities from the THINGS M81 data. An 
approximate sensitivity limit is indicated. HI 
colum ns may become op tically thick at « 5 x 10^^ 
cm-2 (| Allen et al.ll2004 ). 



flux. We evaluate Gq as follows: 



Fig. 2. — Locations of our 27 sources overlayed on 
the GALEX FUV image of M81. 



Go — Fpuv [Dgai/ Phi) /-Fq, 



(3) 



N, 



HI 



7.8x10^ 
.5/50 



In 



1 



106G, 



n ySa J 



-1/2' 



(2) 



where 6 /6a is the dust-to-gas ratio scaled to the 
solar neighborhood value. Nhi (in cm~^) is the 
HI column produced by photodissociation. n = 
riHi + 2n//2(cm^'^) is the total hydrogen volume 
density; this gas is mostly atomic on the surface 
of the GMC and mostly molecular deep inside the 
cloud. Go is a dimensionless parameter measuring 
the strength of the incident UV field. It is scaled 
to twice the equivalent Habing flux (at 1500 A), 
which is appropriate for an FUV-opaque PDR il- 
luminated over 27r sr (jHollenbach fc Tielensl[l999t 
Allen et al . 2004). The latter authors also discuss 
how the flux and its scaling at 1500 A relate to 
the 1000 A flux, which is in the primary range of 
energies of dissociation photons. Go is eq uivalent 
to 0.8 5 times the fi eld strength defin ed in iDraine 
(|l978l ) and used bv lSternber3 (|l988l ). It does not 
include the ambient UV background radiation, so 
it is calculated from a background subtracted UV 



where Dgai is the distance of M81 and pni is 
the distance between the UV source and the lo- 
cation of the atomic hydrogen column. Fq = 
2xFh = 3.11 X IQ-^ ergs cm"^ s~^ A~^ at 1500 A 
(|Allen et al.ll2004 IVan Dishoeck k Blac k 1988). 

Equation [2] can be inverted in order to obtain 
the total hydrogen volume density n, using Nhi, 
Go and S/Sq. The resulting equation is: 

IO6G0 



[exp{NHi{S/So)/7.8 x lO^"} - l]"' (4) 

Equation 0] is valid for PDRs with n < 10^ cm^'^. 

The morpholo gy of the regions identified by 
Allen et"al] (jl997 ) is that of HI patches in the form 
of shells (or parts of shells) partially surround- 
ing young star clusters. Although this morphol- 
ogy appears ubiquitous in M81 and other galax- 
ies, independent confirmation of these HI features 
as PDRs would be very desirable. Such a con- 
firmation could be provided by the presence of 
PAH emission. Although no quantitative relation 
between PAH, UV and HI emission at this scale 
is currently available, the detection of PAHs is 
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average within R25 



R/R5 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
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Fig. 6. — M81 filling factor as a function of galac- 
tocentric radius. The maximum HI column per 
kpc bin, converted to Mq pc^^, is divided by the 
ring-averaged values of M81's radial profile. The 
average filling factor within R25 is 0.21. There is 
no data for i?/i?25 < 0.3. 



thou ght to be predominantly linked to UV emis- 
sion (jUchida et allll998t iLi k Drain j|2002f) . 



Fig. 7. — FUV source fluxes in 10 ergs cm 

s-i A-i. 



4. Method 

In this Section, we explain how the total hy- 
drogen volume density in every candidate PDR is 
determined by investigating the HI structures sur- 
rounding the central FUV source. We hypothesize 
that this source (a cluster of young stars) creates a 
photon dominated region, in which molecular hy- 
drogen is converted into atomic hydrogen. While 
the exact morphology of these regions remains un- 
clear, the distribution of HI suggests that large 
scale photodissociation is taking place. The pres- 
ence of PAH emission in these regions reinforces 
this picture. 

We selected 27 sources based on their well- 
defined compact FUV emission, from faint to 
bright. The sources were chosen to have a range 
of galactocentric radii, although we suspect that 
all candidate PDRs should have similar total 
hydrogen densitie s regardless of their location 



(fSmith et al.ll2000h . . t is significant that we found 
no sources of this type within R/R25 ~ 0.3 and 
brighter than our faintest source ofO.lSx 10~^^ 
ergs cm" -2 A-\ There is one source at a 
smaller radius that has a similar FUV flux as our 
faintest source, but it is extended and has a lower 
contrast than our faintest source. The locations 
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Fig. 8. — Incident flux Go on the HI patches 
chosed for study near the FUV sources. 

of these sources are shown in Figure [2] and given 
in Table □ 

Four example sources are shown in Figure [3] 
These sources vary in FUV flux and galactocentric 
distance. Each source shows the general struc- 
ture of a UV source surrounded or bordered by 
HI with PAH emission ne arby. At la rger scales 
this is clear l y visi ble, see Allen et "all ( ,1997) , or 
Smith et al] i|2000f) for similar plots of MlOl. 

The FUV source flux is determined first. We 
correct for galactic foreground extinction, but not 
for internal extinction. We decided against apply- 
ing an internal extinction correction, d espite the 
availa bility of internal extinction data ( Hill et all 
19951 ). It is assumed that the height of MSl's gas 
disk is of the same order as the size o f the PDRs to 
which our method is most sensitive. IVisser (Il980h 
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Table 1 

Locations and FUV fluxes of candidate PDRs. 



Source 


R. a. (2000) 


DEC(2000) 


Radius 


Ffuv 




Aperture 


Cont. start 


Increment 


no. 


Ch 


m s) 


(d m s) 


(kpc) 


(10"^^ ergs cm-^ 




\ ) 




(10^° cm~^) 


(10^° cm 




1 


9 55 


54.016 


69 


00 


33.79 


4, 1 


1 Qft 
1.00 




17 




o.U 


1.080 


2 


9 55 


41.057 


68 


59 


45.82 


A 7 
4. / 


Q TO 




20 




1.0 


1.170 


q 


9 55 


25.268 


69 


08 


15.24 


A Q 


A CiQ 




10 




O.U 


1.210 


^ 


9 55 


53.182 


68 


59 


04.74 




1 1 Kfi 




24 




Q n 
y .u 


1.190 


c 



9 55 


03.784 


69 


09 


04.41 


0.1 


111 




10 




D.O 


1.120 


ft 
D 


9 55 


10.104 


69 


09 


21.54 


D.l 


i.DO 




10 




O.U 


1.095 


-7 
1 


9 54 


57.442 


69 


08 


49.19 


6.3 


K OA 




15 




D.O 


1.170 


8 


9 56 


15.363 


69 


02 


23.58 


6.4 


o.oU 




16 




A K 
4.0 


1.220 


9 


9 54 


48.655 


69 


06 


00.02 


D. / 


o.oy 




13 




y.u 


1.130 


10 


9 55 


28.266 


68 


58 


47.66 


6.7 


1.57 




13 




4.0 


1.120 


11 


9 54 


53.138 


69 


08 


51.02 


6.8 


4.51 




15 




5.0 


1.180 


12 


9 54 


55.272 


69 


09 


28.97 


7.0 


3.22 




13 




4.5 


1.130 


13 


9 54 


48.210 


69 


03 


22.35 


7.8 


2.91 




16 




11.0 


1.130 


14 


9 55 


03.708 


69 


10 


54.55 


7.8 


1.16 




12 




7.7 


1.080 


15 


9 54 


49.476 


69 


02 


56.82 


7.9 


2.13 




17 




9.5 


1.140 


16 


9 56 


09.346 


68 


56 


48.16 


8.2 


0.88 




13 




6.0 


1.175 


17 


9 54 


39.443 


69 


05 


26.60 


8.3 


3.63 




15 




7.5 


1.150 


18 


9 54 


39.949 


69 


05 


03 44 


8.3 


3.36 




12 




13.0 


1.115 


19 


9 56 


17.179 


69 


05 


42.48 


8.5 


1.87 




13 




20.0 


1.080 


20 


9 56 


17.148 


69 


06 


04.85 


8.8 


4.94 




17 




14.0 


1.135 
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14.667 


68 


57 
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1.87 




13 




15.0 


1.075 


22 


9 55 


25.114 


69 


13 


09.70 


10.9 


2.85 




13 
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1.150 


23 


9 55 


06.623 


69 


14 


08.97 


11.4 


0.66 
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6.0 


1.130 


24 


9 56 


25.503 


68 


53 


43.61 


11.8 


0.53 




15 




9.5 


1.095 


25 


9 55 


08.312 


68 


56 


17.30 


11.9 


0.15 
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9 55 


44.892 


68 


51 


54.80 


14.1 


1.26 




12 
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Fig. 9. 


The distances 


Phi from the UV sources 


Fig. 10 


Derived total hydrogen volume densi- 


to any surrounding patches of HI. 






ties of GMCs in M81 


(n = UHi 


+ 2n//,(cm-3)) 



use a gas disk half height of up to 450 pc, which is 
comparable to the size of our candidate PDRs. In 
that case the extinction from the HI patch to the 
UV source is the same as the extinction of the UV 
source in MSl's disk towards the line of sight of 
the observer. Under these circumstances, the cor- 
responding extinction corrections to the UV flux 
cancel each other out. Using the internal extinc- 



tion correction is therefore unnecessary. 

The radius at which the UV emission blends 
into the background radiation level determines 
the size of the UV source on the image and the 
place to measure the background. The total FUV 
flux within this radius is calculated and the back- 
ground flux is subtracted. Typically, the net flux 
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III 
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Fig. 11. — Various observational limits are shown 
in this plot of Go against n. The numerals are 
keyed to observational selection effects discussed 
in the text. 



is 50% of the total FUV flux. For comparison, 
we calculate the contrast of the source above the 
background level, which we define as the incident 
UV flux on the HI patches G (the unsealed Go) 
divided by the incident background UV flux Gbg- 
The latter is defined as the UV background sur- 
face brightness at the edge of the aperture of the 
UV source, integrated over half the sky (a flux 
on a semi-inflnite slab). Only the net flux of the 
source is used to calculate Go, since the model 
does not include calculating HI produced by the 
ambient UV background. The measured flux is in 
ergs cm~^ s~^ A~^. 

The distance pni from the UV source position 
to any surrounding patches (concentrations) of HI 
is measured next, along with the column density 
of each such patch. A patch of HI is identified as 
a local maximum in column density. See Figure [31 
middle panels. This proce dure differs slightly from 



that of Smith et al 



( 2000l ). who took the first local 
maximum of the ring averaged HI column densi- 
ties around UV sources as a probe of the true HI 
column density. In fact, any local bright spot in 
atomic hydrogen is an equally valid probe of a po- 
tential PDR. They are all likely to be direct views 
onto the surfaces of molecular clouds surrounding 
the UV source. In that case, the HI local maxima 
are a measure of the true local column density 
and each measurement represents (part of) an HI 
cloud. For this reason we decided to use more than 
one measurement of the HI surrounding the UV 



source. We subtract a general background level 
of HI emission that we assume has either not been 
created by photodissociation, or that has been cre- 
ated by the general UV background radiation field 
acting on all GMCs in the area. 

In addition to identifying individual HI concen- 
trations, the maximum measured column densi- 
ties per 1 kpc galactocentric bin are recorded and 
compared to the ring-averaged values from M81's 
radial profile. Dividing the maximum column den- 
sities per bin of individual HI clouds, converted to 
Mq pc~^, by the average surface density per 1 kpc 
bin from M81's radial profile yields a filling factor. 
In the framework of the derived PDR model, this 
filling factor is an estimate of the extent to which 
GMCs fill M81's disk. 

In order to determine pni, we corrected for 
MSl's inclination and position angle assuming 
MSl's disk to be intrinsically circular (see Figure 
2^). From this, the incident flux Go on each sur- 
rounding HI patch is computed using Equation [3l 
The relevant dust-to-gas ratio at the galactocen- 
tric radius of each UV source is taken from Figure 
[TJ The extinction between M81 and the observer is 



taken to be 1.37 mag at 1 500 A, after lAllen et al 



(|l997t) and ISeatonI (|l979l ). The measured FUV 



flux is increased accordingly when calculating Go . 
The incident flux Go, the HI column density and 
the dust-to-gas ratio S/do flnally yields a total hy- 
drogen volume density for each patch according 
to Equation 21 Note that this probe for molecu- 
lar gas does not yield any information about the 
geometry of the gas clouds, but deep inside the 
CMC all available hydrogen is likely to be in the 
form of molecular hydrogen, so n = Sn^f^ ■ 

A PAH map was obtained by subtracting a 
scaled and smoothed version of the 3.6 /zm M81 
im age from the 8 /xm M81 image, as described 
bv IWillner et al. ( 2004 ). This procedure assumes 
that the 3.6 fim image contains purely stellar emis- 
sion, while the 8 /im image contains a mix of stellar 
and PAH emission. We determined the required 
scale factor empirically, by comparing four regions 
of M81 encompassing most of our candidate PDRs. 
Correlation maps of individual pixels were used 
to determine the slope of the appropriate corre- 
lated component. Our four scale fa ctors are con- 
sisten t with the single one used by IWillner et al 



(|2004f ). being within a range of 0.7 and 1. Seven 
sources were not included in these four regions; 
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these sources have comparatively large galactocen- 
tric radii and any 3.6 contribution is negligible. 
In these cases we therefore used the raw 8 /im im- 
age. The relevance of the detected PAH emission 
around the candidate PDR was tested by checking 
for the presence of PAH emission near the location 
of the UV source and determining the separation 
between PAH emission peaks and HI patches. 

Summarizing: determining Nhi, Phi, Fpuv 
and 5/5i) yields n, derived using the picture of pho- 
todissociation. Additionally, the presence of PAH 
emission reinforces the photodissocation assump- 
tion. 

5. Results 

In this Section, we present total hydrogen vol- 
ume densities of the 27 candidate PDRs. Table 
[1] lists the FUV fluxes of the central (cluster of) 
sources. Table [2] lists the corresponding HI col- 
umn density measurements, the resulting incident 
fluxes and finally the total hydrogen volume den- 
sities. First, we'll discuss the HI measurements, 
followed by the FUV results. These results, com- 
bined with the appropriate dust-to-gas ratio from 
Equation [TJ are then used to calculate n. At the 
end of this Section, we discuss the spatial correla- 
tion between HI and PAHs in M81. 

5.1. Atomic Hydrogen 

Each region can have one or more HI measure- 
ments, as was explained in the previous Section. 
The resulting HI column densities and their dis- 
tance to the central FUV source are listed in Table 

H 

Figure [5] shows the HI column densities of all HI 
patches from Table [2] as a function of (normalized) 
galactocentric radius. The general background 
column density level is of the order of 1 x 10^° 
cm~^ and has been subtracted from all measured 
values, for reasons mentioned above (Section [5]) . 
Very little HI radio emission is observed within a 
radius of 3 kpc of the center of M81. The HI col- 
umn densities rise with increasing galactocentric 
radius from 0.3 to 0.6 x i?25 (3.8 kpc to 7.5 kpc), 
although there is significant scatter. 

While the general morphology of our candi- 
date PDRs is the same (FUV emission surrounded 
by atomic hydrogen), the actual geometry of the 
GMCs can not be discerned at the available res- 



olution. The HI structures are barely resolved, 
which makes it practically impossible to deal with 
projection effects properly. This introduces inac- 
curacies in the determination of p^i (Figure [51 
discussed below) and Nhi- 

We can not include HI patches closer than 
about 20 pc (roughly the size of one pixel), or 
right on top of a UV source, even though they 
could have a significant separation in the vertical 
direction. The model can not take a pni of 0. We 
had to reject only one HI patch because it was too 
close to the UV source. The upper limit of pni 
(about 600 pc) is somewhat arbitrary and depen- 
dent on the maximum separation that could still 
conceivably have HI connected to a certain UV 
source. In combination with the lowest measured 
FUV flux, the resuhing Go would be 0.01. 

We corrected for M81's inclination and posi- 
tion angle (Figure H^), assuming M81's disk to 
be perfectly circular in shape. This means we 
used ellipses to measure the deprojected separa- 
tion between the FUV source and the associated 
HI patches. However, pui is still sensitive to pro- 
jection effects, as illustrated in FigureHlD. It shows 
how the measured separation could further be de- 
projected if we knew the position of all HI patches 
in the z-direction of the disk of M81. Therefore, 
the measured separation pni is close to, but not 
exactly equal to the real separation. This will lead 
to underestimating the separation between the UV 
source and its associated HI patches. This in turn 
will mean overestimating Gq and therefore overes- 
timating n as well. At the same time, due to our 
finite resolution, there is a lower limit to the sep- 
aration we can observe. The lowest value of pui 
we measured was 35 pc. 
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Table 2 

hi measurements, incident fluxes and pah correlation. 



Source 
no. 


Phi 
(pc) 


Nhi 
(10^1 cm~^) 


Go 


G/Gbg 


n 


Error 

(%) 


PAH-HI 
< 6 apart 


la 


70 


0.96 


4.18 


3.89 


169 


50 


n 


b 


122 


0.86 


1.37 


1.27 


67 


44 


n 


2a 


209 


2.91 


2.93 


1.41 


5 


88 


y 


b 


297 


2.88 


1.46 


0.70 


9 


87 


y 


3 


175 


2.77 


1.98 


0.47 


8 


85 


y 


4 


52 


4.83 


62.12 


15.20 


27 


140 


y 


5a 


209 


1.78 


0.37 


0.17 


7 


66 


n 


b 


297 


1.81 


0.19 


0.08 


4 


66 


y 


6a 


122 


0.72 


1.62 


0.29 


163 


48 


y 


b 


227 


0.57 


0.47 


0.08 


64 


44 


n 


c 


279 


1.18 


0.31 


0.05 


14 


54 


y 


7 


52 


2.88 


31.93 


8.70 


192 


96 


y 


oa 


1 




1 fiQ 

i.Dy 


U.DO 




79 
1 z 


y 


b 


279 


2.10 


0.66 


0.25 


10 


73 


n 


c 


454 


2.63 


0.25 


0.09 


2 


85 


y 


9a 


52 


2.70 


18.23 


6.58 


151 


93 


y 


b 


122 


2.56 


3.35 


1.21 


32 


85 


y 


10a 


122 


0.89 


1.55 


1.94 


128 


53 


y 


b 


157 


1.07 


0.94 


1.18 


58 


55 


y 


c 


244 


1.11 


0.39 


0.49 


23 


55 


n 


11 


384 


2.13 


0.45 


0.15 


7 


75 


y 


12a 


140 


1.28 


2.43 


0.99 


117 


60 


y 


b 


314 


1.28 


0.48 


0.19 


23 


59 


y 


13 


192 


3.02 


1.16 


0.46 


10 


98 


n 


14a 


192 


1.18 


0.47 


0.48 


31 


61 


n 


b 


244 


1.53 


0.29 


0.30 


12 


67 


n 


c 


349 


1.07 


0.14 


0.15 


11 


59 


y 


15a 


227 


3.27 


0.61 


0.25 


4 


104 


y 


b 


279 


3.09 


0.40 


0.17 


3 


99 


y 


c 


332 


3.45 


0.29 


0.12 


2 


107 


y 


16a 


105 


2.52 


1.18 


0.92 


19 


89 


y 


b 


244 


2.17 


0.22 


0.17 


5 


81 


n 


17 


35 


2.52 


43.88 


16.23 


737 


102 


y 


18 


87 


3.48 


6.51 


1.71 


46 


111 


y 


19 


262 


3.73 


0.40 


0.10 


2 


116 


y 


20 


87 


4.58 


9.56 


2.39 


32 


138 


y 


21a 


105 


2.60 


2.52 


1.36 


59 


97 


y 


b 


367 


2.74 


0.21 


0.11 


4 


98 


n 


22 


506 


3.20 


0.16 


0.05 


3 


113 


y 


23a 


105 


1.57 


0.89 


1.44 


83 


87 


n 


b 


332 


1.11 


0.09 


0.14 


14 


78 


n 


c 


367 


1.81 


0.07 


0.12 


5 


90 


y 
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Additionally, the HI column that is needed for 
Equation|4]is defined as seen along the line of sight 
from the UV source to the HI layer (Figure Ht). 
The observer sees the HI column from a different 
angle. Since the geometry of the HI layer is un- 
known, we are assuming that the observed column 
is nevertheless a reasonable measure of the proper 
HI column. 

The range of observable HI columns is indi- 
cated in Figure [S] With a signal-to-noise ratio of 

« 4 — 5, only regions with column density above 
1 X 10^^ cm~^ can be reliably distinguished as HI 
clouds. {Nhi reaches this level in the region sur- 
rounding source no. 11.) This limits the maximum 
observable total hydrogen volume density. 

At a column density of about 5 x 10^^ cm^^, 
the HI columns may become optically thick 
( Allen et al. 2004) . In our sample, two sources 
have HI columns approaching this value. One of 
these sources has a relatively high UV flux. Nei- 
ther source has an unusual total hydrogen volume 
density. 

Our method is biased towards low density, large 
scale HI regions. The HI from compact high den- 
sity PDRs like the Orion Nebula in our Galaxy 
can not be currently detected at the distance of 
M81 because of the limited resolution of the HI 
map. The high density regions will be of small 
scale size generally and difficult to discern from the 
generally-lumpy HI background. Beam smoothing 
will also reduce the measured value of Nhi in this 
case. Detectable quantities of HI can be produced 
through a combination of a high Go and high n or 
a low Go and low n (recall Equation [21). The high- 
est Go we can detect is limited by the maximum 
observed FUV flux and the minimum measured 
Phi (Equation [3]) . Taking the respective maxi- 
mum and minimum values from our results yields 
a maximum Go of 140. The actual maximum Go 
in our dataset is 62. Obtaining a combination of 
a high Go and a low n is limited by the maximum 
HI column density that is observed. 

Next, we used the HI column densities to calcu- 
late the area filling factor of M81. It is worth men- 
tioning here that Figure [5] wea kly refl e cts th e same 
highs and lows that feature in iBraun ( 1997t ) (their 
Figure 6a) when they used ring averaged values. 
The individual HI patches, presumably indicating 
GMCs, follow the global trend in HI column densi- 
ties. However, the HI columns of individual GMCs 



sh ow too i nuch scatter to reproduce the two peaks 
in iBrauri (1997)'s results. From the individual HI 
column measurements and the global radial dis- 
tribution of HI in M81, the filling factor in Figure 
ini was derived. Table [3] lists the maximum mea- 
sured HI columns per 1 kpc interval, the equiv- 
alent surface densities, the corresponding surface 
densities of MSl's radial profile, and the resulting 
area filling factor. The radial profile has a tilted 
ring correction, appropriate for an inclined disk. 
At a local scale this correction is not appropriate, 
since the morphology of the individual HI columns 
is unknown. Note that we have no sources between 
0.95 and 1.11 x R25 (resp. 12 and 14 kpc). The 
average filling factor within R25 is 0.21. 

5.2. Far Ultraviolet Flux 

Table [1] lists the locations and FUV fluxes of 
the central sources of our 27 selected regions. The 
candidate PDRs were selected to have a range of 
FUV fluxes and galactocentric radii. We did not 
find candidate PDRs in the inner region of M81 
(i?/i?25 < 0.3) with an FUV flux above 0.15x 
10~^^ ergs cm~^ s~^ A~^. Qualitatively, nei- 
ther detectable FUV emission nor detectable HI 
columns are present there. 

Fpuv as shown in Figure [7] shows a weak corre- 
lation with galactocentric radius. No values of the 
FUV flux of individual sources larger than 11.56 x 
10""'^^ ergs cm^^ s~^ were found. At larger 
galactocentric radius, the UV sources become less 
bright. This could be caused by fewer dissociating 
stars per cluster or a change in the stellar pop- 
ulation (an increasing amount of B stars and a 
decreasing amount of O stars). At the lower end 
of the FUV flux values, the minimum measured 
FUV flux values (and Go, Figure [S]) seem to de- 
crease with increasing galactocentric radius. This 
is potentially a confusion effect, if the brighter 
sources raise the ambient UV flux, outshining the 
weaker sources. However, it could also be the ef- 
fect that is responsible for the decreasing maxi- 
mum UV fluxes. At this point we are unable to 
make that distinction. 

Table [5] shows the ratio of the FUV flux to 
the FUV background level [G/Gbg)- A value of 1 
means that the FUV source appears to be equally 
bright as the FUV background as "seen" by the 
HI patch. Values range from a source flux of 1% 
of the background level to over 23 times the back- 
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Table 2 — Continued 



o OLircc 


PHI 


H I 


^0 






Error 


PAH-HI 


no. 


(pc) 


(10^1 cm-") 






(cm-^) 


(%) 


< 6 apart 


24a 


157 


1.81 


0.32 


1.02 


25 


92 


n 


b 


192 


2.12 


0.21 


0.68 


13 


97 


n 


c 


332 


2.44 


0.07 


0.23 


3 


103 


n 


d 


332 


2.12 


0.07 


0.23 


4 


97 


n 


c 


436 


2.28 


0.04 


0.13 


2 


100 


n 


f 


541 


2.04 


0.03 


0.09 


2 


96 


n 


25a 


87 


1.65 


0.29 


0.50 


28 


92 


y 


b 


506 


1.73 


0.01 


0.01 


1 


91 


n 


26a 


192 


1.28 


0.50 


0.87 


107 


100 


n 


b 


297 


1.85 


0.21 


0.36 


27 


109 


n 


27a 


35 


2.51 


9.65 


23.39 


828 


134 


y 


b 


279 


2.83 


0.15 


0.37 


11 


129 


n 


c 


297 


2.20 


0.13 


0.32 


14 


118 


n 


d 


559 


2.75 


0.04 


0.09 


3 


128 


n 



ground level. 

Figure [8] shows the incident flux on all of our 
candidate PDRs. The values of Go show a ten- 
dency to decrease with increasing galactocentric 
radius. Go varies over a range of 4 dex, while the 
FUV flux only varies over 1 dex and shows only a 
very weak trend (Figure [7]). The same holds for 
the values of pni as plotted in Figure [21 ranging 
from 35 pc to almost 600 pc. Taken together, they 
are responsible for the trend in the values of Go- 

5.3. Total Hydrogen Volume Density 

Figure [TO] shows the resu lting total hydro- 
gen v olume density using the IStauffer fc Bothun 
(|l984l ) dust-to-gas ratios (see Equation [1]) . Val- 
ues range between 1 and 200 per cm^, without 
a clea r gradient. Using the iGarnett fc Shield^ 
(|l987l ) dust-to-gas model (not plotted), the re- 
sulting volume densities are in the same range, 
but with larger values in the outer regions of M81 
and smaller values in the inner parts. Even then 
a clear gradient is absent. On average the relative 
error of the resulting n in our method is 88%; see 
Appendix [X] for a more complete error analysis. 

Figure [Tl] illustrates the various observational 
selection effects by showing the range of Gq and n 
that are observable with our method. The differ- 
ent limits are indicated by the Roman numerals I 
through V. These limits arise as follows: 

I. HI column density (sensitivity) lower limit 
of 1 X 10^° cm~^ at a characteristic S/6o of 
1. 

II. HI upper limit of 5 x 10^^ cm~^ at a char- 



acteristic 6/5o of 0.9, due to the HI column 
becoming optically thick. 

HI. The lowest total hydrogen volume density 
that we can observe. It is determined by the 
lowest observable column density divided 
by the diameter of the radio beam and is 
roughly 0.5 cm~'^. This limit is therefore a 
combination of a sensitivity limit and a res- 
olution limit. 

IV. The minimum Gq (0.01) we can reasonably 
use, as mentioned previously. 

V. The maximum Go of 140 that we can obtain, 
also mentioned earlier. 

At the top end we have a theoretical limit of 
roughly 1 x 10^ cm~'^, a combination of the highest 
Go and the lowest HI column density (not shown) . 
These values will not be observed due to points I 
and V above. Moreover, recall that Equation [4] is 
valid only for values of n < 10* cm~'^. More com- 
mon values of Go combined with low HI columns 
yield values of up to 1000 cm~'^. In short, the 
range of values in our results is bounded by ob- 
servational selection, and the apparent trend in 
Figure [11] must be considered a result of that se- 
lection. 

5.4. Polycyclic Aromatic Hydrocarbons 

We now turn to the additional evidence sup- 
porting the idea that the regions we selected are 
indeed PDRs. Morphologically, all sources are 
similar in that HI and PAHs can be found in the 
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neighborhood of the UV emission. This is appar- 
ent from the detailed plots of the candidate PDRs 
in Figure [3l This apparent spatial connection at 
high resolution between the atomic hydrogen and 
UV emission in M 81 was previously reported in 
Allen et al. ( 19971 ) . Almost all of our sources show 
PAH emission nearby. At the largest galactocen- 
tric radii, the PAH emission becomes too faint to 
detect. Table [2] (last column) shows whether the 
measured HI patches have PAH emission peaks 
within 6", or 100 pc. 31 out of 56 HI patches have 
PAH emission peaks within that range. 

Our subtracted PAH map is contaminated by 
non-stellar non-PAH emission (hot dust). In re- 
gions of high radiation intensity, one will observe 
non-stellar emission from HII regions in the 3.6 
//m image, but at a separation of a few arcseconds 
the general view is that the non-stellar emission 
is coming from PAHs (Willner, private commu- 
nication). Despite the contaminations in the re- 
sulting image, we are only interested in a qualita- 
tive comparison, where the presence of emission is 
more important than the actual amounts that oc- 
cur. PAHs are thought to be a tracer of recent star 



vary. 

The cloud densities we find are lower than 
the ra nge of values (30 - 1000 cm^'^) found in 



formation ( Peeters et al. 20041 ). In that case, the 



presence of PAH emission is a further confirmation 
that the physical conditions typical of PDRs are 
indeed present. 

With the exception of 4 sources at higher galac- 
tocentric radius that did not show detectable 
PAH emission (in general, PAH emission levels 
drop with increasing galactocentric radius), all 
UV sources have accompanying PAH emission. 
Owing to differences in excitation conditions, we 
did not expect PAH emission to appear at the 
exact location of HI emission. The majority of HI 
patches, however, has PAH emission associated to 
it, suggesting a close relation between the two. 
A closer investigation of this correlation appears 
desirable, but is beyond the scope of this paper. 

6. Discussion and Conclusions 

The candidate PDRs in M81, which were se- 
lected on the basis of their FUV emission, seem to 
fit the photodissociation model well. Our results 
show no systematically different properties of the 
parent GMCs in different parts of M81. The total 
hydrogen volume density is roughly constant, even 
as the underlying HI, FUV and dust-to-gas ratio 



MlOl (|Smith et al.ll2000f ). The observed values 



of Nhi of individual PDRs are similar to those 
seen in MlOl. The observed HI columns in both 
galaxies abruptly increase at the same normal- 
ized radius (0.3) and appear to decline some- 
what beyo nd 0.7, assuming an R25 of 3 0.3 kpc 
for MlOl (IVila-Costas fc Edmunds! Il992h . The 
range of observed HI columns is also the same. 
The downward tren d in Gg also is consistent with 
Smith et all (|2000l ). when no internal extinction 
correction is applied to their data. FigurefTOlshows 
that the densities of the GMCs in M81 do not ap- 
pear to change with galactocentric radius, consis- 
tent with the MlOl results. 



Kaufman et al. (II999I) modeled the expected 



CO intensity from a PDR for a range of incident 
FUV fiuxes and cloud densities. The range of our 
results would be consistent with modeled CO in- 
tensities below 5 X 10^* ergs cm~^ s~^ sr~^, or 
1 X 10~* ergs cm~^ s^^ sr~^ for the v ast majority 
of sou rces (6.4 K km s^^). Figure 1 in I Allen et al.l 
(|2004l ) shows that for the range of values in Figure 
im the modeled CO intensities are independent of 
n. The low volume densities we find are consis- 
tent with a lack of CO emission in M81 as dis- 
cussed by iKnapen et al.l (J2OO6) , and which those 
authors attribute to insufficient excitation. We 
note that the volume number densities of colliding 
H2 molecules in the GMCs of M81 is 1/2 of the 
values for n calculated here. The mean value we 
have found for n in the GMCs of M81 translates 
into a mean number density for n^f^ ~ lOcm""^, 
well below the values required for coUisional ex- 
citation of CO molecules. In this case, the CO 
emission in M81 is in general subthermal. Their re- 
ported upper limit for Ico is 1-03 K km s~^ for the 
regions they investigated, near our source no. 18. 
The Go and n that we find there are consistent to 
that val ue of Icq (and som ewhat higher) per Fig- 
ure 1 in I Allen et al.l (|2004l ). when beam dilution is 
taken into account. The low CO emission in M81 



is also explored in ICasasola et al.l (|in presa ). who 
again point to a lack of excitation of the molecular 
gas. They find no molecular gas in the nucleus of 
M81. The absence of FUV sources to excite the 
gas is consistent with that finding. 

After accounting for observational and projec- 
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Table 3 

hi surface density per radius interval and resulting filling factor. 



Radius 


Max. HI column 


Equivalent Y^hi 


Radial Xlj/j 


Filling Factor 


(kpc) 


10=1 cm-= 


(Mq pc-=) 


(M0 pc- = ) 




4-5 


2.91 


22.17 


3.88 


0.17 


5-6 


4.83 


36.79 


4.80 


0.13 


6-7 


2.88 


21.94 


4.62 


0.21 


7-8 


3.45 


26.28 


6.41 


0.24 


8-9 


4.58 


34.89 


5.56 


0.16 


9-10 


2.74 


20.87 


5.31 


0.25 


10-11 


3.20 


24.37 


5.42 


0.22 


11-12 


2.44 


18.59 


4.97 


0.27 


14-15 


2.83 


21.56 


2.67 


0.12 



tion effects (see the previous Section), we note 
that in the nearby galaxies in general on which 
our method is applicable, it is most sensitive to a 
combination of low Go and low n. 
Summarizing, our conclusions are: 

• We selected a number of discrete FUV 

sources in M81, which wc consider to be po- 
tential PDRs on the surfaces of the parent 
GMCs. 

• The total hydrogen volume densities of 

GMCs close to clusters of young stars in 
M81 are in the range of 1 < n < 200 cm~^ 
with a geometric mean of 17 cm~^. This is 

approximately ten times lower than GMCs 
in MIDI studied with the same method. 

• The low CMC volume densities are consis- 
tent with a lack of CO emission in M81. 

• MSl's GMCs have a filling factor of « 0.21 
within R25. 

• No candidate PDRs are found in M81 within 
R < 0.3ii25- 

• We have provided a thorough analysis of the 
observational selection effects on our results 
and conclude that, while such effects are 
(necessarily) present in our results, our main 
conclusions as to the range and values of the 
total volume densities of GMCs in M81 are 
not affected. 

• The presence of PAH emission in the neigh- 
borhood of our candidate PDRs lends sup- 
port to our view that the HI patches near 



FUV sources are indeed produced by pho- 
todissociation. PAH emission occurs near al- 
most all UV sources. PAH and HI emission 
coincide in more than half of our sources. 
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A. Error analysis 

To estimate the extent to wliicli tlie different observables influence the result, we performed a basic error 
analysis. It takes into account the precision of the determination of the UV flux F, the Hl-to-UV distance 
p, the dust-to gas ratio 6/ 60 and the HI column densities (N). Note that in this case we subsitute Go with 
its observables F and p. 

Formula |4] for n can be written as 

n^C^ (S/So)-'^^ [exp {N6/S0) - 1]"' , (Al) 
P 

where C is a constant that includes the extinction, the distance to M81, the scaling factor for the flux and 
the other constant factors that went into the model. Also note that N = Nhi/ {7.8 x 10^" cm~^). 



OF J " \dp 
dn \ ^ 9 / dn ^ ^ 



dS J ° \dN 
the standard formula for propagation of uncertainties, where 



(A2) 



^ = - (A3) 
dF F ^ ' 

dn n , . , , 

^ = -2- A4 
dp p 

n [exp {N5/5o) - 1]"^ A^exp {NS/Sq) (A5) 
-n[exp {N5/5o) - 1]~^ • 

(5/Joexp(iV5/'5o) (A6) 



dn 
dN 



The square root of the result, (T„, is the standard deviation, or the error. The dust-to-gas ratio has 
a disproportionate influence on the error, since it occurs twice in the equation and is also the parameter 
that is the least well known. The relative precision of the measurements of F and N is roughly 1%. The 
separation has an uncertainty in the order of a half a pixel, or about 9 pc. Finally, the scatter in the fit to 
the dust-to-gas ratio is about 0.2. 

Taken together, these figures introduce relative errors to the result of about 88% of the measured value 
of n on average, or approximately a factor of 2. This means that for example n — 1 ± 0.8 if ^ = 80%. 
Negative values of n are unphysical. The spread in our resulting values of n seems to refiect this and we 
expect the accuracy of our results to be similar. 



16 



REFERENCES 

Allen, R. J., 2004, "Photodissociation and the 
Morphology of HI in Galaxies" , in Penetrating 
Bars Through Masks of Cosmic Dust: The Hub- 
ble Tuning Fork Strikes a New Note, eds. D.L. 
Block, I. Pucrari, K.C. Freeman, R. Groess, 
& E.K. Block (Springer; Dordrecht), 731-748 
( |astro-ph/0407004| . 

Allen, R. J., Atherton, P. D., Tilanus, R. P. J., 
1986, Nature, 319, 296 

Allen, R. J., Knapen, J. H., Bohlin, R., Stecher, 
T. P., 1997, ApJ, 487, 171 

Allen, R. J., Heaton, H. I., Kaufman, M. J., 2004, 
ApJ, 608, 314 

Braun, R., 1997, ApJ, 484, 637 

Casasola, V., Combes, F., Bettoni, D., Galletta, 
G., accepted for publication in A&A 

van Dishoeck, E. F., Black, J. H., 1988, ApJ, 334, 
771 

Draine, B. T., 1978, ApJS, 36, 595 

Freedman, W. L., Hughes, S. M., and 13 coau- 
thors, 1994, ApJ, 427, 628 

Garnett, D. R., Shields, G. A., 1987, ApJ, 317, 82 

HiU, J. K., et al., 1995, ApJ, 438, 181 

HoUenbach, D. J., Takahashi, T., & Tielens, 
A. G. G. M. 1991, ApJ, 377, 192 

Hollenbach, D. J., Tielens, A. G. G. M., 1999, Rev. 
Mod. Phys., 71, 173 

Issa, M. R., MacLaren, I., Wolfendale, A. W., 
1990, A&A, 236, 237 

Kaufman, M. J., Wolfire, M. G., Hollenbach, D. J., 
Luhman, M. L., 1999, ApJ, 527, 795 

Kcnnicut, R. C., et al., 2003, PASP, 115, 928 

Knapen, J. H., Allen, R. J., Heaton, H. I., Kuno, 
N., Nakai, N., 2006, A&A, 455, 897 

Li, A., Draine, B. T., 2002, ApJ, 572, 232 

Mihalas, D., Binney, J., 1981, Galactic Astron- 
omy: Structure and Kinematics 2nd ed. (Free- 
man: San Francisco) 



Peeters, E., Spoon, H. W. W., Tielens, 

A. G. G. M., 2004, ApJ, 613, 986 

Seaton, M. J., 1979, MNRAS, 187, 73 

Smith, D. A., Allen, R. J., Bohlin, R. C., Nichol- 
son, N., Stecher, T. P., 2000, ApJ, 538, 608 

Stauffer, J. R., Bothun, G. D., 1984, AJ, 89, 1702 

Sternberg, A., 1988, ApJ, 332, 400 

Uchida, K. I., Sellgren, K., Werener, M., 1998, 
ApJ, 493, L109 

de Vaucouleurs, G., de Vaucouleurs, A., Cor- 
win, H. G., Buta, R. J., Paturel, G., Fouque, 
P., 1991, Third Reference Catalogue of Bright 
Galaxies (New York: Springer) 

Vila-Costas, M. B., Edmunds, M. G., 1992, MN- 
RAS, 259, 121 

Visser, H. C. D., 1980, A&A, 88, 149 

Walter, F., Brinks, E., de Blok, W. J. G., Thorn- 
ley, M. D., Kennicutt, R. C, 2005, ASP Conf. 
Proc, 331, 269 

Willner, S. P., Ashby, M. L. N., Barmby, P., Fazio, 
G. G., Pahre, M., Smith, H. A., Kennicutt, 
R. C, Jr., Calzetti, D., Dale, D. A., Draine, 

B. T, and 8 coauthors, 2004, ApJS, 154, 222 

Wolfire, M. G., Tielens, A. G. G. M., & Hollen- 
bach, D. 1990, ApJ, 358, 116 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.2. 



17 



